Critical free energy and Casimir forces in rectangular geometries 
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We study the critical behavior of the free energy and the thermodynamic Casimir force in a 
LjJ -1 x L block geometry in 2 < d < 4 dimensions with aspect ratio p — L/L\\ on the basis of the 

O(n) symmetric if 4 lattice model with periodic boundary conditions and with isotropic short-range 
interactions. Exact results are derived in the large - n limit describing the geometric crossover from 
film (p = 0) over cubic (p = 1) to cylindrical (p = oo) geometries. For n = 1, three perturbation 
approaches in the minimal renormalization scheme at fixed d are presented that cover both the 
central finite-size regime near T c for 1/4 < p < 3 and the region well above and below T c . At bulk 
T c we predict the critical Casimir force in the vertical (L) direction to be negative (attractive) for 
a slab (p < 1), positive (repulsive) for a rod (p > 1), and zero for a cube (p = 1). Our results for 
finite-size scaling functions agree well with Monte Carlo data for the three-dimensional Ising model 
by Hasenbusch for p = 1 and by Vasilyev et al. for p = 1/6 above, at, and below T c . 

PACS numbers: 05.70. Jk, 64.60.-i, 75.40.-s 



I. Introduction and overview 



In the theory of finite-size effects near phase transitions, 
the study of critical Casimir forces [l[ has remained on 
a highly active level over the past two decades 0] ■ Close 
to criticality and for sufficiently large confining lengths, 
such forces are predicted to exhibit universal features in 
the sense that, for isotropic systems with short-range in- 
teractions, their scaling functions depend only on the ge- 
ometric shape, on the boundary conditions (b.c), and 
on the universality class of the system. For anisotropic 
systems of the same universality class (e.g., lattice sys- 
tems with noncubic symmetry such as anisotropic su- 
perconductors 0L. however, two-scale factor universal- 
ity [i[ is absent [5l-flc| and the critical Casimir forces 
are nonuniversal as they depend on the lattice structure 
and on the microscopic couplings through a matrix of 
nonuniversal anisotropy parameters. This implies that 
the Casimir amplitudes at bulk T c depend, in general, 
on d{d + l)/2 — 1 nonuniversal parameters in d dimen- 
sional anisotropic systems with given shape and given 
b.c. [fHH. This prediction is readily testable by Monte 
Carlo (MC) simulations of anisotropic spin models; it was 
also noted in [7[ that experiments in anisotropic super- 
conducting films 0] could, in principle, demonstrate the 



nonuniversality of the critical Casimir force . A corre 



sponding nonuniversality of the Binder cumulant [1, [I?} 
has been predicted 0, 0] and has been confirmed by MC 
simulations for the anisotropic Ising model 

In our present paper, the focus is on isotropic systems 
with periodic b.c. for the (n = 1) Ising universality class 
on the basis of the 0(n) symmetric ip 4 lattice model. 
Theoretical studies (beyond mean-field theory) of the 
Casimir force in such systems have been restricted to 
oo d ~ 1 x L (film) geometry within the e = 4 — d expansion 
above and at T c [l9l - [2lT |. The most interesting region, 
however, is the region below T c where the scaling function 
of the Casimir force displays a characteristic minimum as 
detected by MC simulations [TU, [22j • A basic difficulty 
in treating an infinite film system (for n — 1) below T c 



is the existence of a film transition at a separate critical 
temperature T c jn m < T c just in the region of the mini- 
mum of the Casimir force (23j . A quantitative theory of 
the corresponding dimensional crossover between differ- 
ent critical behavior of the d dimensional bulk transition 
at T c and the d — 1 dimensional film transition at T c jn m 
constitutes an as yet unsolved problem, except for the 
case of the Gaussian model Q . 

In the present paper we circumvent the problem of di- 
mensional crossover by studying a finite iy _1 x L block 
geometry with finite aspect ratio p = L/L\\ [24]. This 
geometry includes slab (0 < p < 1), cubic (p = 1), and 
rod-like (p > 1) geometries (Fig. 1). The practical rele- 





FIG. 1: (Color online) Three-dimensional Ly x L block ge- 
ometries with aspect ratio p = L/L«: (a) slab (p < 1), (b) 
rod (p > 1). Arrows: critical Casimir force (|2.15[) . Our the- 
ory (in Sees. Ill - V) predicts that, for isotropic systems with 
periodic b.c. and short-range interactions, Fcasimir,s at bulk 
Tc is negative (attractive) for a slab, positive (repulsive) for 
a rod, and zero for a cube (p — 1). 

vance of the slab geometry is based on the facts (i) that 
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all experiments and all MC simulations have necessarily 
been performed at finite p rather than at p = 0, (ii) that 
the shape of the finite-size scaling function of the Casimir 
force depends on the aspect ratio p only weakly in the 
regime p <C 1 [lj| [22j , and (iii) that the singularity of the 
free energy and the Casimir force at T c ju m for p = is 
only very weak such that recent MC data [l6j could not 
detect this singularity. (By contrast, the logarithmic di- 
vergence of the specific heat for p = near T c ju m should 
be well detectable.) This justifies to describe the main 
features of the film system above and below T c , to a good 
approximation, by a finite slab geometry with small but 
finite aspect ratio. As an interesting by-product of our 
theory, the dependence on the aspect ratio for larger p 
is obtained which permits us to describe the geometric 
crossover from slab (p <C 1) over cubic (p = 1) to rod- 
like (p ^> 1) geometries. The geometric crossover from 
block to cylindrical (p — oo) geometries has been studied 
earlier near first-order transitions by Privman and Fisher 
[25j . In this context we note that systems with p 3> 1 are 
of experimental interest in the area of finite-size effects 
near the superfluid transition of 4 He [26|, [27j. Further- 
more, finite-size theories for block geometries are directly 
testable by MC simulations. 

The finite-block system is conceptually simpler than the 
infinite film system because of a considerable technical 
advantage: For finite < p < oo, the system has a dis- 
crete mode spectrum with only one single lowest mode, in 
contrast to the more complicated situation of a lowest- 
mode continuum in film (p = 0) or cylinder (p = oo) 
geometry. This opens up the opportunity of building 
upon the advances that have been achieved in the de- 
scription of finite-size effects in systems that are finite in 
all directions [71 |28H31| . It is not clear a priori, however, 
in what range of p such a theory is reliable since, ulti- 
mately, for sufficiently small p < 1 or sufficiently large 
p > 1, the concept of separating a single lowest mode 
must break down. Therefore, as a crucial part of our 
theory, it is necessary to provide quantitative evidence 
for the expected range of applicability of our theory at 
finite p. This is one of the reasons why we also consider 
(in Sec. Ill) the large - n limit whose exact results turn 
out to be helpful in estimating the range of validity of 
our approximate results for n = 1. 

As we shall present three different perturbation ap- 
proaches with different ranges of applicability for the case 
ii = 1 we give here a brief overview of our strategy. The 
basic physical quantity is the singular part of the excess 
free energy density f^ x . On the basis of previous work 
7, 32] it is expected that, for our isotropic system with a 
finite volume £j| x L, there exist three different regimes 
(a), (b), and (c) for the finite-size critical behavior of the 
excess free energy f^ x . These different regimes corre- 
spond to the three regions in Fig. 2 that are separated 
by the dashed lines. 

(a) The regime well above T c that includes an expo- 
nential size dependence ff* ~ exp(— L/£ e +) or ff* ~ 
exp(-L||/£ e +) for large L/^ e+ > I and L||/£e+ > 1 
with £ e+ being the exponential ("true") bulk correlation 
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FIG. 2: Schematic plot of the asymptotic part of the L" 1 ^ —t 
plane, with t = (T — T c )/T c , for the (p 4 model with n = 1 and 
with isotropic interactions on a simple-cubic lattice with lat- 
tice spacing 5 in a L^T 1 x L block geometry with periodic 
boundary conditions and finite aspect ratio p = L/Ln. In the 
central finite-size region (above the dashed lines), the lowest 
mode must be separated whereas outside this region (below 
the dashed lines), ordinary perturbation theory is applica- 
ble. Above the shaded region, universal finite-size scaling is 
valid for isotropic systems. In the large - L regime at t 
(shaded region), the exponential form of the size dependence 
violates both finite-size scaling and universality because of a 
non-negligible dependence on the lattice spacing 5. An anal- 
ogous plot applies to the L7, 1 — t plane. 



length @,IH,|33| above T c ; in this regime, f* x is expected 
to tend to zero in the high-temperature limit at finite L 
and L|| or in the limit of large L and L\\ at fixed temper- 
ature T > T c . 

(b) The central finite-size regime near T c that includes 
the power-law behavior ff ~ L~ d or ff° ~ LJT for large 

L at fixed L/£± ,0 < L/£± < 0(1) and for large L\\ at 
fixed Lii/£± ,0 < Lu/£± < 0(1) above, at and below T c 
where £± is the second-moment bulk correlation length. 

(c) The regime well below T c that includes an exponen- 
tial size dependence ~ exp(— L/£ e _) or ~ exp(— Lii/£ e _) 
for large L/£ e - ^ 1 an d L||/£ e - ^ 1 with £ e _ b eing the 
exponential (" true" ) bulk correlation length 0, 1331 |3~H 
below T c ; in this regime, ff is expected to have an ex- 
ponential decay towards a finite value —V~ x In 2 in the 
low-temperature limit at finite volume V [32| and to tend 
to zero in the limit of large volume at fixed temperature 
T <T C . 

For a description of the cases (a) and (c) , ordinary pertur- 
bation theory with respect to the four-point coupling uq 
of the ip 4 theory is sufficient. This ordinary perturbation 
approach is applicable to the regions below the dashed 
lines in Fig. 2. This approach will be presented in Sec. 
VI. For the case (b), a separation of the lowest mode and 
a perturbation treatment of the higher modes is neces- 
sary 0, [28l - [3l| . This approach is applicable to the region 
between the dashed lines in Fig. 2 which we refer to as 
the central finite-size regime. In Sees. IV and V we treat 
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the case (b) in 2 < d < 4 dimensions on the basis of the 
<p 4 lattice model in the minimal renormalization scheme 
at fixed dimension d (35[. We consider a simple-cubic 
lattice with isotropic short-range interactions. We shall 
demonstrate that our different perturbation approaches 
complement each other and that the corresponding re- 
sults match reasonably well at intermediate values of the 
scaling variables. As will be shown in Sees. V and VI, 
these results agree well with Monte Carlo data for the 
three-dimensional Ising model by Hasenbusch for p = 1 
(36j and by Vasilyev et al. for p = 1/6 [16] above, at, and 
below T c . 

We shall see that in all regimes (a)-(c), universal finite- 
size scaling [37j of isotropic systems is valid, except for 
the regions that are indicated by the shaded areas in Fig. 
2 where the exponential form of the size dependence vio- 
lates both finite-size scaling and universality because of a 
non-negligible dependence on the lattice spacing a. The 
boundaries of the shaded areas are not sharply defined; 
they are approximately determined by L ~ 24^ /a 2 . 
This issue will be discussed in Sec. VI D. 



is the dimensionless partition function. The bulk free 
energy density per component divided by k B T is obtained 
by 



f b (t) = lim lim f{t,L,L\\). 

L—too L II — >oo 



(2.6) 



The film free energy density per component divided by 
ksT is obtained by taking the limit Lu — ¥ oo at fixed 
finite L (i.e., p — > 0) 



ff limit, L) 



lim 



f(t,L,L\\), 



(2.7) 



corresponding to an oo d 1 x L geometry. Our model also 
includes the limit L — > oo at fixed finite L\\ (i.e., p -4 oo) 
corresponding to an Lf~ 1 x oo geometry which we shall 
refer to as cylinder geometry, 



fcyl(t,L\\ 



lim f(t,L,L«). 

L—>oo 



(2.8) 



The excess free energy density per component divided by 
k B T is 



r(t,L,L«)=f(t,L,L«)-f b (t). 



(2.9) 



II. Model and basic definitions 



We start from the O(n) symmetric </? 4 lattice Hamiltonian 
divided by fc^T 



For the finite LfJ 1 xL system we define the Casimir force 

per unit area and per component in the d th (vertical) 
direction (Fig.l) as 



H 



N N K 

£(^ 2 + U0 (^ 2 ) + £^r^-^ 



(2-1) 



r (T) = r 0c + aot, t=(T-T c )/T c , 



(2.2) 



with ao > 0, uo > where T c is the bulk critical tem- 
perature. The variables (ft = v?( x i) are ^-component 
vectors which represent the internal degrees of freedom 
of N particles on N lattice points Xj of a d-dimensional 
simple-cubic lattice with lattice constant a and with pe- 
riodic boundary conditions. The components ipf^ , [i = 
1,2,..., n of tfii vary in the continuous range — oo < 



-d-l 



x L 



<Pi ' < oo. We consider a finite rectangular Ljj 
block geometry of volume V = V^ 1 L with the aspect 
ratio 



L 



(2.3) 



This block geometry includes the shape of a cube (p = 1), 
of a finite slab (0 < p < 1), and of a finite rod (1 < p < 
oo) (Fig.l). 

The free energy per component and per unit volume di- 
vided by ksT is 

f^L,^) = -(nV)- 1 In Z(t,L,Lu) (2.4) 



where 



Z(t, L, in 



11 jjn(2-d)/2 



cxp(-iJ) (2.5) 



Casimir 



[t,L,Ln 



d[L.r(t,L,L {l )} 
dL 



(2.10) 



where the derivative is taken at fixed Lu. There exist, 
of course, also Casimir forces in the d — 1 horizontal di- 
rections. Our approach is well suitable to calculate such 
forces. This will not be performed in this paper. 

An important simplification of our model Hamiltonian 
(|2.ip is the assumption of a rigid lattice with a rigid 
rectangular shape representing an idealized model system 
with a vanishing compressibility. The same assumption 
is made in lattice models on which previous Monte Carlo 
simulations of the Casimir force are based (see [1EE3). 
As a consequence, the number N of particles and the 
length L are directly related hy N = L?r l L/a d . Thus 
the derivative with respect to L (at fixed Lu, fixed a 
and at fixed couplings Kij and u Q ) in (|2.10[) is equiva- 
lent to a derivative with respect to the number of hor- 
izontal layers of the lattice, i. e., number of particles. 
Such a definition of the Casimir force is appropriate when 
the ordering degrees of freedom (particles in fluid films 
EH EBl or Cooper pairs in superconducting films Q) 
can move in and out of the film system without signif- 
icantly changing the mean interparticle distance in the 
system. The definition (|2.10p is, however, not appropri- 
ate for systems with a fixed number of ordering degrees 
of freedom. It appears that this is reason why it was 
claimed in [l2j that the Casimir force "is not a measur- 
able quantity for magnets". For similar claims see [38| . 
There exist, however, (long-ranged) critical fluctuations 
of elastic degrees of freedom coupled to the order param- 
eter in condensed matter systems with a finite compress- 
ibility (such as magnetic materials [lj], alloys and 
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solids with structural phase transitions [T3|) which give 
rise to L dependent critical Casimir forces that are not 
contained in a description based on rigid-lattice models 
of the type (|2.1[) . A description of such thermodynamic 
Casimir forces is provided by coupling the variables tpi 
in (|2.ip to the elastic degrees of freedom [l3|, EE > m 
which case the free energy density fit, L, L\\,N) depends 
on both the length L and the number N of particles as 
independent thermodynamic variables. Such a model is 
relevant for the calculation of an L dependent elastic re- 
sponse to the critical Casimir force (e.g., an L dependent 
contribution to magnetostriction). In the case of a com- 
pressible system whose number N of particles is fixed, 
the L dependent thermodynamic Casimir force is given 
by 



(t,L,L h N) 



d[Lr-(t,L,L lb N)} 
dL 



(2.11) 



where now the derivative is taken at fixed t,L\\, and N. 
In such systems, anisotropy effects on the critical Casimir 
force are expected to play an important role. Eq. (|2.11j) 
complements our arguments presented in [Tl| . As we con- 
sider the thermodynamic Casimir effect as a finite-size ef- 
fect, our definition (|2.1ip does not include the bulk part 
of the total thermodynamic force — d[Lf(t, L, L\\ , N))/dL 
which may give rise to measurable elastic bulk effects 
(such as a bulk contribution, e.g., to magnetostriction). 
Here we shall not further discuss this extension to sys- 
tems with a finite compressibility which is beyond the 
scope of the present paper whose focus is on isotropic 
incompressible systems. 



For small \t\, the bulk free energy density (|2.6I) can be 
uniquely decomposed into singular and nonsingular parts 



fb(t) = /&,,(*) + h,ns{t). 



(2.12) 



For large L/a, L\\/a and small \t\, a corresponding as- 
sumption is made 4] for 

f{t, L, Ly) = /,(*, L, L„) + f na (t, L, L|, ) (2.13) 

and for the excess free energy 

f ex (t,L,L l{ ) = /r it, L, + fZ it, L, L||) (2.14) 

where / ns (t, L, Ly) and f^{t, L, Lii) are regular functions 
of t and where f ns (t,L,L») remains regular in the bulk 
limit, / ns (t,£,L||) ->■ fb,ns(t), whereas f s (t,L,L\\) -> 
fb,s{t) becomes singular in this limit. It has been as- 
sumed [37| that, for periodic b.c, f ns (t, L,L\\) is inde- 
pendent of L and Ly, i. e., that it is identical with the 
regular bulk part fb,nsif) of fb(t). We do not know of a 
general proof of this property; it appears to be valid for 
the if 4 theory in the presence of short-range interactions 
Kij but not in the presence of long-range correlations 
0. The critical behavior of Fcasimir can be calculated 
from its singular part 



Fcasimir,si^i L, Ltt) 



0[LCT,£,£|l 
dL 



(2.15) 



For the <p 4 lattice model with the interaction given in 
(|2.19[) below it is expected that f%%it,L,Ln) = 0, thus 



Fcasimir = Fcasimir.s which is consistent with our results 
in Sees. Ill - VI (see also the remark after Eq. (4.36) of 
0)- 

Our main goal is to study the case n — 1 at fixed finite 
aspect ratio < p < oo including extrapolations to the 
film and cylinder geometries. For comparison and as a 
guide for our extrapolations we shall also consider the 
exactly solvable limit n — > oo. 

For periodic boundary conditions, the Fourier represen- 
tations arc 



,ik-x., 



and 



K i j = K{x i -* j ) = N-^> 



(2.16) 



K(k) , (2.17) 



where the summations ^ k run over N discrete vectors 
k = (k\, k2, . . . , kd) with Cartesian components k a = 
2irm a /L\\, a = 1, 2, • • • , d — 1, and kd — 2Trm,j/ L, mp = 
0, ±1, ±2, • • • , in the range —ir/a < kp < n/a, f3 = 
l,...,d. In terms of the Fourier components the Hamil- 
tonian reads 

1, 



^ 2 



H = V- 1 > - [r + Jif (k)]^(k)^(-k)) 



u V-* $30(k)0(p)][0(q)0(-k p q)] (2.18) 

kpq 

where ^^(k) = 2[A"(0) — -f^(k)]. We assume a simple 
ferromagnetic nearest-neighbor interaction 



5K(k) = — [1 - cos(afe Q )] 



which has the isotropic long-wavelength form 
6K(k) = k 2 + Oik 4 ) . 



(2.19) 



(2.20) 



Thus it is appropriate to define a single second-moment 
bulk correlation lengthf ± above (+) and below(— ) T c 
(see, e.g., Eq. (3.4) of [7[). As a reference length that is 
finite for both n = 1 and n — oo we shall use the asymp- 
totic amplitude £o+ of the second-moment bulk correla- 
tion length above T c 



(2.21) 



For small \t\, the asymptotic bulk power law is /b, s (i) = 
At \t\ du . Due to two-scale factor universality for isotropic 
systems Q, this can be written as 



/MW-i i A-/A+)Q 1 i^\tn d 



for T >T C , 
for T <T C , 

(2.22) 

with a universal constant Q\ and the universal ratio of 
the specific-heat amplitudes A~ /A + . 



The finite-size scaling form of the singular part of the free 
energy density is, for isotropic systems in the asymptotic 
critical region \t\ < 1, L/a > l,LJa » 1, 



f s it,L,L«) = L- d Fii,p) 



(2.23) 



where F(x, p) is a dimensionless scaling function with the 
scaling variable 



x = t(L/Co+) 1/v - 



(2.24) 



The bulk part F*(x) of F{x,p) is obtained from (jUg) 
and (|2.23p in the limit of large \x\ as L d f s (t,L,L\\) -)• 
F^(x) with 



hi**" 



for T > T c , 



F *( x )- < ! (A-/^+)Qi7i for T < if (2 ' 25) 



This implies the scaling form 



^(i,p) = F(i,p)-F b ± (i) 



(2.26) 



(2.27) 



Together with (|2. 15[) , this leads to the scaling form of the 
critical Casimir force for systems with isotropic interac- 
tions 



Fcasimir,s{t,L,L\\) — L X(x,p) 

with the scaling function 



(2.28) 



v ox op 

(2.29) 

These scaling functions have finite limits for p — >• at 
fixed L and at fixed i corresponding to film geometry, 



ffilm,s(t,L) — L d Ffu m (x), 



Ffum{x) = lim F(x,p), 

p— >0 



(2.30) 



(2.31) 



f™(t,L,L l{ ) = L^ d *<*(x h p), (2.37) 

F C asimir,s(t, L , L \\) =^ d E(x\\, P ), (2.38) 

with the scaling functions 

$(x h p) = p- d F(x ll p 1 /»,p), (2.39) 

$™(x h p)=p- d F e *(x ll p 1 /»,p), (2.40) 

E(x h p) = p- d X(x^/",p) 

= -^ x {x hP ) + (1/p) ^ PJ {2A1) 

It turns out that they have finite limits for p —> oo at fixed 
L|| and fixed xii corresponding to cylinder geometry, 

fcyl,.(t> L \\) = L i d *c V l(x\\), (2-42) 



$ cyi (i||) = lim $(x {l ,p), (2.43) 



*% l {x\\)=]im* e *(x h p), (2.44) 

" p— >oo 



S cyZ (J,|) = lim S(i||,p) = -*^(i||). (2.45) 

Quantitative results for the various scaling functions will 
be presented in Sec. Ill for n — oo and in Sees. V and 
VI for n = 1. 



F% m (x) = ]imF™(x,p), 
1 p— >o 



Xfum(x) = lim X(x,p) 

p— s-0 

3 v ox 



(2.32) 



(2.33) 



Note that v denotes the bulk critical exponent and x is 
the scaling variable with respect to bulk criticality, not 
with respect to film criticality. 

In a rod-shaped geometry with p ^> 1, a representation 
of the scaling form in terms of the length L« and of the 
scaling variable 



i,l =<(i||/eo + ) 1/ ^, 



(2.34) 



rather than in terms of x, is more appropriate. Because 
of 



x = Xup 1 ^ 11 



(2.35) 



we obtain from ([235) , (|2~25|) . (j^lj) . (g^BJ) , and (|2~27|) 

/.(*, £, = , p), (2.36) 



We recall that all finite-size scaling forms given in this 
and the subsequent sections are not valid in a small part 
of the asymptotic region for large L (or large L«) at fixed 

t ^ in the L 1 /" — t plane (or lV v — t plane) above and 
below T c (corresponding to the shaded regions in Fig. 
2) where exponential nonscaling terms exist that depend 
explicitly on the lattice constant a. This issue will be 
discussed in Sec. VI. 



III. Large - n limit 

A. Free energy density 

Generalizing Eq. (6.30) of Ref. to Lm 1 x L geometry 
we obtain for the free energy density per component in 
the limit n — > oo at fixed uqu 



f(t,L,L«) = lim [-(nV)~' L kiZ(t,L,Ln 



ln(27r) 
2h d 



"'" x J,J ■.i ? 5>Kx- 1 +«£(k)])a a ] ( 3.i) 



16uo n 
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where Z(t, L, Ly) is defined by (|2.5I) and \(t, L, Lh 
determined implicitly by 



In the film limit p — > at finite L, we obtain pK(p 2 z) 
(^) 1/2 , 



X 1 = ro 



Auqh 
V 



^ix-' + SKik)}- 1 . 



The condition \ 1 = for bulk (L — > oo,Ln 
criticality yields the critical value of tq as 



roc 



i n J 



7T I a 



4u n / Jif(k)- 1 , 



n 



2tt 



(3.2) 

> oo) 

(3.3) 



Ad 
2s 



j,\2 



zP(£) 



P(J) 



I 



J j, film 



(P 2 ) = (4tt 2 )- j / d«2 i_1 eacp 



— 7r/a 



{(tt/z) 1 ^-^)} 



£o,/iim(P(£) ), 

(3.10) 

P 2 Z\ /7T\(d-l)/2 

(3-11) 



Eqs. d3H]) and (02]) are valid for T >T C (r > r 0c ) and 
T < T c (r < r 0c ). For T > T c , the quantity x(t,L,L||) 
represents the susceptibility per component. 



B. Exact scaling functions 



In the following we present exact results for the finite- 
size scaling functions F, $, F ex , <& ex , X, and S in 2 < 
d < 4 dimensions. We rewrite (|3.ip and p. 21) in terms of 
?"o — J"0c = oo^ aud assume large L/a, large Lu/a, small 
\r -r c\/a 2 < 1 and |r -r 0c |L 2 < 0(1), \r Q -r 0c \Lj < 
0(1). Evaluation of the sums in p. II) and (|3.2[) (see App. 
A) leads to the scaling form of f s (t, L,L\\), (|2.23[) and 
(f2T24|) . where v = (d - 2)" 1 and 



where P(x) is determined implicitly by 



P d ~ 2 = x - — GutimiP 2 ). (3.12) 
Ad 

For n — > oo, a finite film-critical temperature < 
T c jiim(L) < T c exists only in d > 3 dimensions whereas 
Tcju m (L) = in d < 3 dimensions. Eqs. (j3~TU]) - (|3~T2]) 
are valid in the asymptotic region near bulk T c . 



In the cylinder limit p 
K(p 2 z)^l,P(x,p)/p^P cyl (x ll ) 



oo at finite we obtain 



Ad 



5||P cy ;(i||) 2 - -j F ^(*||) c 



1 



:0O,cy/(Pcyz(J||) 2 ), 



(3.13) 



4u nA d \ 1 /(rf-2) 



V eao I 



with the geometric factor 



A d = 



T(3 - d/2) 
2<*-2 7r d/2(d_2)" 



(3.4) 



(3.5) 



4tt2 



For an arbitrary finite shape factor < p < oo we find 
the finite-size scaling function 





x {Wz)W-[K{z)]^}, 
where P cy i(xu) is determined implicitly by 

x 



cyl 



pd-2 



'II 



F(x,p) 



2e 



iP(x,p) 2 - jP(i,p) d 



+ \Go(P(x,p) 2 ,p)- 



e 

T d 



Gl,cyl(Pcyl) 



(3.14) 
(3.15) 



(3.6) 



OO 



d-l 



x (tt/z)^ 2 - P K(p 2 z) K(z) (3.7) 



(3.8) 



where P(x,p) is determined implicitly by 



P 



d-2 



£-^-Gi(P 2 , P ) 

Ad 



with 



K(z) = exp(-zm 2 ). 



(3.9) 



The earlier result of Eqs. (6.32)-(6.34) of for cubic 
geometry is obtained from p.6[) - (|3.8[) by setting p = 1 . 



In this cylinder limit, the system has an infinite exten- 
sion only in the d th direction, i.e., it is essentially one- 
dimensional, thus no finite critical temperature exists in 
the cylinder case at finite L\\. Eqs. (|3.13p - p,15[) are 
valid in the asymptotic region near bulk T c . 

The bulk part Fp(x) of F(x,p) in the large - n limit 
is given by (|235l) with Qi = (d - 2)A d /[2d(4 - d)] and 
A~ /A + = 0. Correspondingly, the bulk part of $(x||,p) 
is $^(i||) = P fc ± (iy). From (|2~27|) and (gHQjl we then 
obtain the scaling functions F ea: and Q ex of the excess 
free energy density. 

The scaling functions X and 5 of the Casimir force are 
obtained from F ex and $> ex according to ([2T25]) . (j2351) . 
(I2.41j) . and (|2.45[) . All scaling functions are shown in 
Figs. 3 and 4 for several values of p in three dimensions, 
illustrating the crossover from film geometry (p = 0, dot- 
ted lines in Fig. 3) over cubic geometry (double-dot- 
dashed lines) to cylinder geometry (1/p — 0, dotted lines 
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-0.2 



-0.4 ~ 



-0.6 



(a) 



p = 0, film 

p = 1/4 

p = 1/2 

■-■ p = 3/4 
p = 1, cube 






-0.2 



0.4 
ex 

-0.6 



(a) 



0.8 


_ s' 


1/p = 0, cylinder 




t s ' 


1/p = 1/4 


-1 


_ #»• ' 


1/p = 1/2 




■ — ■ ■ 1/p = 1, cube 



(b) 



l/p= 


0, cylinder 


— l/p = 


1/4 


--- l/p = 


1/2 


l/p = 


1 , cube 




FIG. 3: (Color online) Scaling f unctio ns (a) F ex (x,p), p^T) . 
[pTrjf - IpO]), and (b) X(x,p), QTZty , as a function of x = 
t{L/^o+) 1 ^ v in the large-n limit in three dimensions for film 
geometry (p = 0, dotted lines), for slab geometry (solid lines: 
p = 1/4 , dashed lines: p = 1/2, dot-dashed lines: p = 3/4), 
and for cubic geometry (p = 1, double-dot-dashed lines). For 
x —5- — oo, the curves in (a) with p > diverge logarithmically 
towards — oo whereas the p — curve in (a) has a finite low- 
temperature limit -0.191. All curves in (b) have finite limits 
for i -> -oo as given by (|3.2ip and (|3.23p . 



x hP ), magi, 

l/u 



FIG. 4: (Color online) Scaling functions (a) $ 
and (b) S(x«,p), (|2.41|l . as a function of = t(.L||/£o+ 
in the large-n limit in three dimensions for cylinder geometry 
(1/p = 0, dotted lines), for rod geometry (solid lines: 1/p = 
1/4, dashed lines: 1/p = 1/2), and for cubic geometry (1/p = 
1, double-dot-dashed lines). For arii — > — oo, the curves in (a) 
with 1/p > diverge logarithmically towards — oo whereas 
the 1/p = curve in (a) has a finite low-temperature limit 
-0.719. All curves in (b) have finite limits for x« — > — oo as 
given by ([3722]) and lpOi)l . 



in Fig. 4). We see that, for 0{— 10) < x < oo and 
0{— 10) < x\\ < oo, the scaling functions for slab (p < 1) 
and rod (1/p < 1) geometries, respectively, provide a 
reasonable approximation for the scaling functions (i) for 
film geometry if the shape factor p is sufficiently small, 
(ii) for cylinder geometry if the inverse shape factor 1/p 
is sufficiently small. This is not true, however, in the low- 
temperature limit x — > — oo and jcn — > — oo, respectively 
(see the following subsection). 



C. Monotonicity properties at fixed p 

For fixed p, F ex (x,p) and & ex (x\\,p) are monotonically 
increasing functions of x and x\\, respectively [see Figs. 
3(a) and 4(a)]. They vanish exponentially fast for x — > oo 
and a5y — >• oo and have logarithmic divergencies towards 
—oo for x — > — oo and xu — > — oo, respectively, for finite 



< p < oo. 

To derive the latter property, consider the quantity P 
as determined by (|3.8[) . It is finite and positive for 
— oo < x < oo and vanishes for x — > — oo. More specif- 
ically, the function Q\ has the divergent small-P 2 be- 
havior ^i(P 2 ,p) ~ -p d ~ x P- 2 [see (|A~8l) in App. A]. 
According to (|3.8|) . this implies that P 2 vanishes as 
P 2 w eA' d 1 p d ~ 1 (— x)~ x for x — > — oo. Thus the behavior 
of F(x, p), (I3.6p . for large negative x is given by 

F(x, p) = F°°(x, p) « -\? d - X + \Go(P(x, p)\ p) 

(3.16) 

for —x 1. The function Qq has a divergent small-P 2 
behavior as given by (|A.4I) in App. A. The resulting 
logarithmic divergency is 

P) « ^ (^P) - \p d ^ + \C { P ) 

(3.17) 
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with Co(p) given by (IA.6I) . For finite < p < oo, Eqs. 
(|2.35|) and (I2.40[) imply a corresponding logarithmic di- 
vergency of <& ex (£\\, p) for afii — >■ — oo. 



By contrast, we shall find a nonmonotonic dependence 
of the scaling functions F ex and <I> ea: on x and acii for 
the n — 1 universality class for finite < p < oo in 
the central finite-size scaling regime described in Sec. V 
below. 

For the film system in the large-n limit, we confine our- 
selves to the case d = 3. We find from ()3.1ip that 
Ql,film(P 2 ) ~ (47r)- 1 lnP 2 for small P 2 and that P 2 



vanishes as P oc e x for x 



-oo. This implies that 



Fex 
Him 



(x) has a finite value in the low-temperature limit 
[compare Fig. 3(a)] 



lim F]f lm {x) = - Go,fUm(0) = -0.191 for d = 3. 



(3.18) 



In the cylinder system in the large-n limit, we find 
from (|3~T4l) that gi, cy i(P 2 yl ) « -ciP~] for small P cyi 
with a = {Air)- 1 / 2 J °° dyy-^e-v > and P cyi « 
Ci£j47 1 (— in)" 1 for sen — > — oo. This implies that, for 
2 < d < 4, nas a /™^ e value in the low- 

temperature limit 



lim 



(3.19) 



with | 0o )C! /i(°) = -0-719 for d = 3 [compare Fig. 4(a)]. 

In contrast to F ex and $ ex , the scaling functions X and 
S turn out to be nonmonotonic functions of their scal- 
ing variables x and xii, respectively, in an intermediate 
range of p where p ~ O(l) [see Figs. 3(b) and 4(b)]. 
In this range, X exhibits a change of sign near x = 0: 
The Casimir force changes from an repulsive force be- 
low T c to an attractive force above T c . Especially for 
p = l,d = 3, this change of sign occurs exactly at T c 
where X(0, 1) = 0, S(0, 1) = 0. In the range p < 1/2, 
X < is a monotonically increasing function of x. In 
the range /9>3/2,S>0isa monotonically decreasing 



function of x 



Above T c , X and S have an exponential decay towards 
zero as functions of x 3> 1 and x\\ 3> 1, respectively, 
as follows from the exponential decay of F ex and Q ex . 
Below T c , the scaling functions X and 5 have finite values 
in the low-temperature limits x — > — oo and scm —J- — oo, 
respectively, for all — oo < p < oo, unlike the divergent 
behavior of F ex and & ex for finite p. To derive the low- 
temperature limit of X we use (|2.27p and (|3.6[) to rewrite 
(|2~2^|) as 



X(x,p) = F±(x) + — 



e 



—xP — 



1 



(d-i)g (P 2 ,p)- P 



d :p ° 
dg (P 2 ,p) 



dp 



(3.20) 



with P{x,p) 2 determined by (f378|) . For P 2 -)• 0, the 
divergent parts of the last two terms cancel each other 



which leads to a finite limit 
1 



lim X(x,p) = --p 



d-l 



x^ — oo 



l-\im\(d~l)G (P 2 ,p) 

dg (P 2 , P y 



d P 



(3.21) 



with a nontrivial p dependence. Similarly we obtain a 
finite limit 



lim E(x h p)=p- d lim X(x,p). (3.22) 

||— V — oo x— > — oo 



This is in contrast to the exponential decay of X and 5 
towards zero for x — > — oo and jcn — > — oo, respectively, 
for the n = 1 universality class that we shall find in Sec. 
VI below. 





-0.5 




-1 




-1.5 


o 






-2 




-2.5 



-3.5, 




n=infinity 

n=l 

• MC Ising model 



0.5 



1 

P 



1.5 




FIG. 5: (Color on line) Critical amplitudes (a) F(0,p), ^2~Tf\ . 
and (b) X(0,p), (|5.10(l . at T c in three dimensions as a func- 
tion of the aspect ratio p in the large-n limit [thin lines, from 
(|3.6[) ] and for n = 1 [thick lines, from l|5.1[l ]). For n = oo, 
F film (0) = F(0,0) = -0.153 and X film (0) = X(0,0) = 
2P/« m (0) = -0.306. At p = 1, X(0, 1) vanishes for both 
n — oo and n = 1. Monte Carlo data of the d — 3 Ising model 
for p — 1 by Mon [4lJ [full circle in (a)] and for p — by 
Vasilyev et al. lj| [triangle in (b)]. See also Fig. 7. 



For film and cylinder geometries in the large - n limit, 
the low-temperature limits are 

lim X fllm (x) = 2FS m (-oo) = Goju m (0) = -0.383 

x— >~ OO 

(3.23) 
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-0.5 



Q- 
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n=infinity 












• MC Ising model 
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1/ P 
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FIG. 6: (Color online ) Critical amplitudes (a) $(0,p) (f2739]) 
and (b) H(0, p), (|5.11|) . at T c in three dimensions as a function 
of the inverse aspect ratio 1/p in the large-n limit [thin lines, 
from (13.6[) ] and for n = 1 [thick lines, from (JET}]). For n = oo, 
$ c „i(0) = $(0,oo) = -0.329 = -H cyi (0) = -H(0,oo). At 
1/p = 1, E(0, 1) vanishes for both n = oo and n — 1. Monte 
Carlo data of the d — 3 Ising model for p = 1 by Mon [4ll ] 
[full circle in (a)]. See also Fig. 8. 



for d = 3 [see Fig. 3(b)], in agreement with the earlier 
result for the spherical model [22|, E(3| , and 

lim E cyl (xn) = -$cj,i(-CX5) = -- Qo ,cyl(0) (3-24) 
xii— > — oo Z 



towards zero for p — > and for 1/p — ► 0, respectively. 
For comparison, corresponding curves are also shown for 
the n = 1 universality class (thick lines in Figs. 5 and 6) 
that will be derived in the subsequent Sections. On the 
basis of these results we are led to our hypothesis that 
the monotonicity properties mentioned above are valid 
not only for n — oo but are general features of the free 
energy and the Casimir force (for periodic b.c.) that are 
valid for all n in the whole range 1 < n < oo. 



IV. Perturbation theory for n = 1 in the central 
finite-size regime 



In this and the subsequent sections we confine ourselves 
to the case of a one-component order parameter. 



A. Perturbation approach for the free energy density 



The basic ingredients of our perturbation approach for 
Ljj -1 x L geometry are similar to those developed pre- 
viously for cubic geometry 0- The starting point is a 
decomposition of the variables tpj = $ + <jj into the low- 
est (homogeneous) mode amplitude $ and higher-mode 
contributions cr,-, 



$ = V~ x tp{Q) = N- 1 



3 



(4.1) 



(4.2) 



Correspondingly, the Hamiltonian H and the partition 
function Z are decomposed as 



H = H + H($,a), 



(4.3) 



H (r Q , u , V , $ 2 ) = V ( ^r $ 2 + w $ 4 ) , (4.4) 



for 2 < d < 4, with -| g cy i(0) = 0.719 for d = 3 [see 
Fig. 4(b)]. 



D. Monotonicity properties at fixed temperature 



a) = d d l [(y + Guo&' Z )v- + 4u $<t? + u a; 



(4.5) 



From Figs. 3 and 4 we also infer monotonicity properties 
at fixed x and xn, respectively, i. e., at fixed tempera- 
ture. For fixed x, F ex (x,p) is monotonically decreasing 
with increasing p whereas X(x,p) is monotonically in- 
creasing with increasing p. For fixed Xii, both $ e:E (x||,p) 
and S e2: (xii , p) are monotonically decreasing with increas- 
ing 1/p. This monotonicity is demonstrated by the thin 
lines in Figs. 5 and 6 at bulk T c (x = 0, xm = 0). These 
lines also exhibit a monotonic change of the curvature 



yV2 



exp{- [# + f($ 2 )]}, (4.6) 



T($ 2 ) = - In 



n I da{k) 



exp 



-H(p,v)\ , 
(4.7) 
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where <r(k) = <^?(k) for k ^ 0. We shall calculate the par- 
tition function and the free energy by first determining 
r($ 2 ) by means of perturbation theory at given $ and 
subsequently performing the integration over $. Since 
exp[— r(<I> 2 )] is proportional to the order-parameter dis- 
tribution function, which is a physical quantity in its 
own right, we shall maintain the exponential form of 
exp[— r($ 2 )] without further expansion. 

The decompositions (|4.3I) - (|4.7p and the perturbative 
treatment of the higher modes are reasonable as long 
as there exists a single lowest mode that is well sepa- 
rated from the higher modes. This is, of course, not the 
case in the film limit p — > and in the cylinder limit 
p — )■ oo where the system has a lowest mode continuum 
and where a revised perturbation approach would be nec- 
essary. In Sec. V below, a quantitative estimate will be 
given to what range of < p < oo our perturbation 
approach is expected to be applicable. 

Since the details of the perturbation approach for 
f(t, L, Lii) are parallel to those presented in [7| for cubic 
geometry we directly turn to the result. Our perturba- 
tion expression for the bare free energy density reads 



/= - 



N - 1 1 

— — ln(27r) + -S (r oh , L, p) 



If 1 

-In ds exp(-- 



— In 

2V 



eff 2 

y s 



V l ' 2 wl ff 



6u Mo5i(r L,i, p) 



36ulMgS 2 (r h,L,p) 



(4.8) 



with 

eff 

Vo 



V^ 2 u~ 1/2 {r [l + mu S 2 (r 0L ,L,p)} 
UuoS^rauL^p) + U4u 2 M 2 S 2 (r 0L , L, p)](4.9) 



The positivity of tql > for all ro permits us to apply the 
theory to the region below T c . For finite V, Mq and r L 
interpolate smoothly between the mean-field bulk limits 
above and below T r 



lim M 2 =M 2 mf 

V— loo J 



and 



lim r L 



for r > , 

-ro/(4uo) for r o < , 

(4.15) 

;° ^ ro ^, (416) 

-2r for r < , v ' 



respectively. In the bulk limit, Eq. (|4.8[) correctly con- 
tains the bare bulk free energy density in one-loop order 
[i.e., up to 0(1)], 

ft = ^ + ~ / ln{[r + ^(k)]a 2 },(4.17) 



fb = ~ r M mf + u M ™f 



1h(2tt) 
2a d 



+ - I ]n{[-2r + 6K(k)]a'} 



(4.18) 



above and below T c , respectively. [For the symbol J 

k 

see ((331) .] In the derivation of (|4.8p . f($ 2 ) has been 
expanded around Mq in powers of <& 2 — Mq up to 
0(($ 2 - Mq) 2 ). Furthermore, an expansion with respect 
to uq at fixed tol has been made and has been truncated 
such that terms of 0(u^ 2 ) are neglected. For a discus- 
sion of the order of the neglected terms see @, Hi], |42| . 



B. Dependence on the aspect ratio p 



<" = «o 1/2 [1 + 18u 5 2 (r 0L , L, p)] . (4.10) 
Here <Si(roL, L, p) denote the sums over the higher modes, 

5 (r 0L , L, p) - In { [r 0L + 5K(k)\ ~a 2 } , (4.11) 



S m (r QL ,L,p) = ^Y,{[ r ^ + S K(k)}} m (4.12) 

k#0 

for m = 1,2. The temperature dependence enters 
through the parameter 

roL(r ,u ,V) = r + l2u a M 2 {r ,u ,V) (4.13) 

as well as through the lowest-mode average 



Eqs. (|4.8|) - (|4.10l) are identical in structure with Eqs. 
(4.26) -(4.28) of [7fl where a cubic geometry was consid- 
ered. The new point of interest here is the dependence 
of the bare free energy density f(t,Ln,L) on the aspect 
ratio p = L/Lu. The p dependence enters (i) through the 
volume 



V = Lfr x L = L d p l - C \ 



y a =r Q (L d p 1 - d /u ) 1 / 2 , 



(4.19) 



(ii) through the lowest-mode average 

M 2 (r , uo, V) = {L d p 1 - d u )- 1 ' 2 tfadfo), (4.20) 



(4.21) 



J d$ $ 2 exp[-H (ro,u ,V)} 

M 2 (r , u , V) = ^ . (4.14) 

J d$ exp[-H (r ,u 0l V)] 



Mvo) 



J ds s 2 exp(— \yos 2 — s 4 ) 
o 

oo 

/ ds exp(-iy s 2 - s 4 ) 
o 



(4.22) 
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and (iii) through the higher-mode sums Si(r ^, L, p), 
(|4.11[) and (I4.12j) . In the regime of large L/a, large Ln/a, 

small < (r 0L ) 1/2 a < 1 and fixed < L(r 0L ) 1/2 < 0(1), 
< L^roh) 1 ^ 2 < 0(1), these sums are evaluated for fi- 
nite 0< / o<ooin2<<i<4 dimensions as (see App. B 
of and our App. A) 

S (r L,L,p)= f\n{[SK(k)]d 2 } + r L /[^(k)]" 1 



2A d (ro L ) rf/2 , 1, ( L 2 \ , 1 - p d - 1 _ 2 



de 



; In , 

L d V a 2 47r 2 



+ -j-jJo(rohL 2 ,p), 



(4.23) 



.s',(/„l.i./m= /[^(k)]- 1 - ^( roh r-^ 



1 - (T\ 2 ~ d 

-^(roL)- 1 + { -^ y I 1 (r ^,p), (4.24) 



5 2 (r 0L ,L,p) = ^ (d-2)(r 0L )- £ / 2 + ^-^(r 0L )- 2 
+ J^2y l2 ^ L ^P) ( 4 - 25 ) 

with 



J (x 2 ,/9)= / dyy 1 exp [-a; 2 2 //(47r 2 )] 



{(7r/y) d / 2 - [p K{p 2 y)] d - 1 K(y) + l} - e -»] ,(4.26) 



/ m (x 2 ,p) = / dy 2/ m - 1 expha:V(47r 2 )] 



x{ [p K(p 2 y)] d - 1 K(y) - (ixjyfl 2 - 1} (4.27) 
for m = 1, 2. [For if (y) see (|3"1?|) .] 

C. Bare perturbation result 



It is appropriate to rewrite the free energy density /, 
|, in terms of rg — rg c where 



r 0c = - 12ii JiSkik)}- 1 



(4.28) 



is the critical value of r up to O(u ). The resulting 
function is denoted as f(jQ — ro c , "o, L, p, Kij, a). As we 
are interested only in the singular part we subtract the 
non-singular bulk part up to linear order in tq — r$ c , 

/£>(r - r 0c , 1^,5) = + \! ln{[8K(k)]~a 2 } 



2a: 



[SKQt)]- 1 . (4.29) 



The remaining function 

Sf(ro -r 0c ,u a ,L,p,Kij,a) = f(r - r 0c ,u , L, p, K id ,a) 

- f£(ro-r c,K itj ,a)(4,30) 

has a finite limit for a — > at fixed r$ — rQ C in 2 < d < 4 
dimensions, 

lim Sf(ro~r 0c ,uo,L,p,Ki d ,d) = 6f{r ~r 0c ,u ,L,p) , 

a— >-0 

(4.31) 

where we have assumed the interaction (|2.19[) . (For the 
justification of taking the limit a —> see the remarks 
after Eq. (4.36) of Q.) The function (|4~3T]) still contains 

(2) 

a non-singular bulk part fn S ' (ro — r$ c , uq) proportional to 
(ro — i'0c) 2 - It is convenient to subtract this non-singular 
bulk part later within the renormalized theory in the 
asymptotic critical region as described in Subsect. E. Our 
perturbation result for the function Sf(ro — ro c , uq, L, p) 
as derived from P~g]) - P~T4]) and (l^7T^|) - (|4"T3T1) reads 



18u^M 4 



Sf(r Q - r 0c ,u ,L,p) = 



A d 

e/2 
r 0L 



PL 

Ad 



(ro - r 0c ) s 
4e 



{OO 
-/-Mn | exp(-iy e// (p)z 2 -z 4 ) 



"5- 



27T^ // (p) 
L e/2 p (d-l)/2 

3u M 2 L 2 



+ 2 J o( r ohL 2 , p) 



2tt 2 

d-1 f ^ 



w T 2 x 9w 2 M 4 L 4 r2 



L d 

r 0L 



In 



4tt 4 
Wo e// (p)r OL i 2 



L e/ 2/9 (d-l)/2 27r 



6m M 2 



OL 



% // (p) = 
1 - /" X 

+ 12u 

+ 144u 2 M 2 
L 



L d/2 (l-d)/2 



1/2 



(ro - r 0c ) 



1 + 18u 



(4.32) 



A d (d-2) 



2e r, 



e/2 
OL 



U 



16tt 4 
-4,, 



1 - /- 1 



£ r, 



(2-d)/2 ' 
OL 

Aj(d-2) 



roL 



4tt 2 



■h(rohL 2 ,p) 



1-P 



d— 1 



2e r, 



s/ 2 

0L 



OL 



16tt 



hirohL ,p) 



(4.33) 



25 



'OL 



+ 



1-p' 



d-1 



'OL 



16tt 4 



hirohL ,p) 



A d (d-2) m - e/2 

(4.34) 
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where now and Mg are abbreviations for 

roh{ro - r Qc , u ,L, p) = r - r 0c + 12m M 2 (4.35) 

and 

M 2 (r - roc, u , L, p) = (LV^«o)" 1/2 Mm), (4.36) 
with 

yo = (ro-r 0c )(L d p 1 - d /u ) 1 / 2 . (4.37) 



Our Eqs. (|4.32[) - (|4.37[) are applicable to some finite 
range of < p < oo and contain Eqs. (4.37) - (4.42) of 
as a special case for p = 1. They are not applicable 
to the film (p —> 0) and cylinder (p -» oo) limits below 
bulk T c . 



D. Minimal renormalization at fixed dimension 



As is well known, the bare perturbation form of Sf re- 
quires additive and multiplicative renormalizations. As 
the ultraviolet behavior of Sf does not depend on the as- 
pect ratio p, the renormalizations are the same as those 
described in 0] in terms of the minimal renormalization 
at fixed dimension 2 < d < 4 [35] ■ The adequacy of 
this method in combination with the geometric factor 
A d , (|3.5p . has been demonstrated in pj for the case of 
cubic geometry. Since the aspect ratio p is not renormal- 
ized we apply the same renormalizations to the present 
bare expression for Sf, (|4.32p . The details are parallel to 
those in Q which justifies to turn directly to the renor- 
malized form of Sf. It is defined as 



f R (r,u,L,p,p) = 5f(Z r r,p, e Z u Z v 2 A d x u,L,p,) 



-a e r 2 A d A{u,e) 



(4.38) 



where r and u are the renormalized counterparts of 
fo — r 0c and ug. For the Z-factors Zi(u,e) and the addi- 
tive renormalization constant A(u, e) we refer to Q. The 
inverse reference length p is chosen as /i -1 = £o+ where 
£o+ is the asymptotic amplitude of the second-moment 
bulk correlation length above T c . 

The critical behavior is expressed in terms of a flow pa- 
rameter l(t,L,p) that is determined implicitly by 



2,2 



r h (l)=p 2 l 



(4.39) 



The reason for this choice of the flow parameter is given 
after (|4.45[) below. The dependence of I on t, L, and p en- 
ters through the function which is the renormalized 
counterpart of tql- It is given by 



r L (Z) = raL (r(l),l e ^Afu(l),L,p) 

= r(l) + 12\(piyA^ 1 u(l)L- d p ( - d - 1) 



1/2 



(4.40) 



with 



y(l) = r(l)p-H- 2 {LpV) d l 2 p^l 2 A X l 2 u{iy x l 2 (4.41) 



where #2(2/) is defined by (I4.22p . The effective renor- 
malized quantities r(l) and u(l) are defined as usual [35| . 
Both and y(l) depend on t, L, and p. The t depen- 
dence originates from r(l) which depends on t through 
its initial value r(l) = r = at with a — Z r (u,e)~ x a§. 

The effective renormalized counterparts of j/q (p) ano - °f 
Wg (p) are given by 

y^f{hp) = {l l iL)^p^-^A 1 J\{l)- 1 ' 2 



IpL, p 




xR 
and 

w eff {Up) 
where 



+lU{lpL)- d ' 2 p( d - i y 2 A- 1,2 u{lfl 2 d 2 {y{l)) 

n,(0 



p?X- 



-,lpL,p) 



(4.42) 



u(iy 1/2 [l + 18u(l)R 2 (^0-, IpL, p 



(4.43) 



Rl(q,p,p) = e- l q[l - q^ 2 } + A^\l - p^q^p^ 
+ pC- 2 (4ir 2 A d )- 1 I 1 (qp 2 ,p), (4.44) 

R 2 {q, P ,p) = -e- 1 [l- q -z' 2 ]- l - q -z/ 2 

+A- 1 (l-p d - 1 )q- 2 p- d + p%16^A d )- l I 2 (qp 2 ,p) , 

(4.45) 



with I m defined by (|4.27p . The dependence of the func- 
tions Ri on the ratio r^(l)/(p 2 l 2 ) is the reason for the 
choice (|4.39p of the flow parameter. It ensures the stan- 
dard choice in the bulk limit both above and below T c 

M 



p 2 l\ = r{l+) 



lim lim n 2 l 2 = i 2 % '+ 
l->oo f.^ I p 2 lt = -2r(l 



for T > T c , 
for T < T c , 

(4.46) 

and implies pi oc L~ x p^ d ~ x ' / d for large finite V at T = T c . 



After integration of the renormalization-group equation 
(see Eqs. (5.6) and (5.7) of [7]), the renormalized free 
energy density attains the structure 



f R (r,u,L,p,n) = f R (r(l),u(l),lp,L,p) 
A d r{l) 2 



+ 



2{lpY 



B(u(l')){ 



exp 



dl" } di 
1 



(4.47) 



where B{u) and (, r {u) are well known field-theoretic func- 
tions of bulk theory 0, |Hj]. From (|4~3"2]) and (|4T55j) we 
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derive the first term on the right-hand side of (|4.47|) as 
fR(r{l),u(l),lfi,L,p) 

= - A d (lp) d /(4d) + 18u(l)L- d p d - 1 k9 2 (y(0)] 2 



{OO 
" / _1 ln / dz ex P [ " \ 



27TA 1 J 2 w^f(l,p) 
{IpLy/ipid- 1 )/ 2 

3(^L) £ /V( rf - 1 )/ 2 u(0 1/2 



9(lnL) e p d ~ l u(l) 
i d I 2 



[^(2/(0)] 2 / 2 av^ 2 ,p) 

[ipLY^p^-^m-K 
^u(i)^(ipL)- d / 2 P ^y 2 A- 1/2 My(i)) 

-36u(l)(lpL)- d p d - 1 A^ 1 [^ 2 (y(/))] 2 ). 



(4.48) 



E. Finite-size scaling function of the free energy 
density 



It is straight forward to show that the asymptotic form 
(|2.23j) of the singular part f s of the free energy density is 
obtained from (|4.47[) , in the limit of small I <C 1 or 

/ — > as 

= (4.49) 



with the scaling variable x, (|2.24|) . In this limit we have 
«(/) -> u(0) = u*, r(l)/(p 2 l 2 ) -> Q* 1 1- 1 /", 

y(l) y = ZQ* {p,lL)- a /^)p^- d )l 2 A l J 2 u*- 1/2 , 

(4.50) 

and plL — >■ Z = l(x,p) where the function l(x,p) is de- 
termined implicitly by 



+ 12tf 2 (y) = p^/ 2 l d ' 2 A X J 2 u*- 1/2 , (4.51) 



y= xQ* p^-^A^u*- 1 / 2 . (4.52) 

These two equations also determine y = y(x,p). In 
(|4.50[) - (|4. 52[) we have used the hyperscaling relation 
2 — a = dv. The factor Q* = Q(l,u*,d) is the fixed 
point value of the amplitude function QCLu, d) of the 
bulk correlation length above T c 0, HH, SI] . Further- 
more we have, in the small - I limit, 



„»// 



(l,p)^ W(x,p) = 



-1/2 



1 + 18 u*R 2 (l,l,p) 



y eff (l,p) Y(x,p) = j^ p H-d)/2 A i/2 u ,-i/2 



x \ Q*x r 1/lJ 



144 



l + l8u*R 2 (l,l,p) + 12u*iii (1,1, p) 



1/2 



**l-d p (d-i) A -i M y )R2{1 J jp) y (4.54) 



For I <C 1, the last integral term in (|4.47[) contains both 
a contribution oc t 2 l~ a / u to the singular finite-size part 
f s (t, L, L») and a contribution oc t 2 to the nonsingular 

bulk part /, ls ; b of Sf mentioned after (|4.31| (sec also the 
comment on Eq. (6.8) of *7]). This nonsingular part will 
be neglected in the following. 



Eqs. g47]), dHH]), and (j4~49l) - (j4~54| lead to the finite- 
size scaling function 



F(x,p) = -A d 



l d vQ* 2 x 2 l- a l" , ^ 
— + — B(u*) 

'2irA]/ 2 W(x,p)\ 



+ 18u*p d - L [My)] -Trln, 
- p A ~ x ln y exp [ - \y{x, p)z 2 - z A ] 

— OO 

2 27r 2 Ay 2 



^d— 1 



9 Z e u*p' 



4^ 4 A d 



[Mv)] 2 h(J 2 , P ) 



+ (l-p d - i ) -In 



^ /2 iy(x,p)/ rf / 2 
27rp(d-i)/2 



-6«* 1/2 r d / 2 ^ 1)/2 ^ 1/2 ^&) 



(4.55) 



with A d , J ,I m and tf 2 defi ned in ([3~5)l . (|4T2d) , (|4T2"7) , 
and (|4.22l) , respectively. Eq. (|4.55[) is the central analytic 
result of the present paper for the case n = 1. It is valid 
for 2 < d < 4 in the central finite-size regime (between 
the dahed lines of Fig. 2), i.e., in the range L 3> a, 
Ln 3> a and < \x\ < 0(1) above, at, and below T c for 
finite p. For the special case p = 1, Eq. (|4.55p is identical 
with Eq. (6.10) of 7]. For d = 3, Eq. (|4~35|) reduces to 
Eq. (9) presented in [2J] . It incorporates the correct bulk 
critical exponents a and v and the complete bulk function 
B(u*) (not only in one- loop order). The Borel resummed 
values of the fixed point value u* fH[, of B(u*) and 
of Q* [3l|, [H, HH in three dimensions are given after Eq. 
(|5.5p below. There is only one adjustable parameter that 
is contained in the nonuniversal bulk amplitude £o+ of 
the scaling variable x. For finite L and Lu, f s (t,L,L\\) 
is an analytic function of t near t = 0, i.e., F(x,p) is an 
analytic function of x near x = at finite /?, in agreement 
with general analyticity requirements. 

The bulk part F b (x) of F(x,p) is obtained from (|4.55j) 
(4.53) in the large - \x\ limit. It is represented by (|2.25p . with 
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the universal bulk amplitude ratios 
Qi = - A d Q* d » 



(4.56) 



A~ 

= 2 



dv l/(64u*) +l/(4d) + 81u*/64 + i/B(«*)/(8a) 



l/(4d) + z/S(u*)/(2a) 



(4.57) 



given by Eqs. (6.19) and (6.20) of 0. We then obtain 
from (|435j) the scaling function F ex (x,p), (|2~27|i . of the 
excess free energy density which determines the scaling 
function X(x, p) of the Casimir force according to (|2.29l) . 
By definition, the functions F ex (x,p) and X(x,p) have 
a weak singularity at x = arising from the subtraction 
of the bulk term F^(x). 



V. Quantitative results in three dimensions in the 
central finite-size regime 
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A. Amplitudes at T c and monotonicity hypothesis 
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Of particular interest is the finite-size amplitude at T c , 
F(0,p)=(l8-^)u*p d - 1 [M0)} 2 

2 n V ^V*- 1 " 2 / 

oo 

- p^ 1 In J dz exp [ - \Y c {p)z 2 - z 4 ] 

— oo 

+ \j^lp)-^h(Jlp)-^- A h(Slp) 



1 



+ (l-p— ) -In 



4 /2 w c ( P )i 



d/2 



2 7 r j o( d - 1 )/ 2 

where l d/2 = \2u* 1/2 p^ 1 ^ 2 A d 1/2 tf 2 (0) and 

W(0, P ) = W c (p) =u*~ 1/2 |i + 18u*JZa(l,i c ,p) 



- S (5.1) 



(5.2) 



F(0, / j)=y c (p) = 144uM 2 (0)|i?i(l,f c ,p)+ J R 2 (l i r c ,p)| 

(5.3) 

withtf 2 (0) = r(3/4)/r(l/4) and 
l 2 - d 



Ri(l,l c ,p) 



A-K 2 A d 



h<p e ,p) + Aj\\-f-i)X-\ (5.4) 



L p) = -\ + iSa 12 ^ p) + Ad ' (1 ~ pd ~ 1]I c d 



(5.5) 



FIG. 7: (Color online) Critical amplitudes (a) F(0,p), ((237)) . 
and (b) X(0, p), ()5. 10p . at T c in three dimensions as a function 
of the aspect ratio p for n = 1 [thick lines, from (|5.ip ]~) and 
in the large-n limit [thin lines, from (|3.6[) ]. The maximum 
-0.1636 of the n = 1 line in (a) is at p max = 0.2470. The 
dashed lines are the extrapolations of the n = 1 lines from 
p ~ Pmax to p = corresponding to film geometry. The 
dotted lines represent (|5.1[) in the regime p < p ma x where our 
perturbation theory is not applicable. MC estimate for the 
d=3 Ising model from [161 ] (triangles), e expansion results for 
n = 1 from 19] (squares) and from [20J (diamonds). See also 
Fig. 5. 



For the application to three dimensions we shall employ 
the following numerical values @, HH, |45| : A 3 = (47r) -1 , 
v = 0.6335, u* = 0.0412, Q* = 0.945, B{u*) = 0.50, and 
a — 2 — iv = 0.0995. At T c in three dimensions, the p 
dependence of the flow parameter is given by 

l c = [12(4 7 ru*) 1 / 2 r(3/4)/r(l/4)] 2 /V /3 = 2.042 p 2 / 3 . 

(5.6) 

~ 2 

The p dependence of the integrals Jq(I c ,p), (|4.26[) . and 
I m (l 2 ,p), (|4.27p . needs to be computed numerically. The 
resulting amplitudes F(0, p), X(0, p), $(0, p), and S(0, p) 
as determined by ([5TT]) . (|5TTU|) . (|2"39"|) , and ([5~TTj) arc 
shown by the thick lines in Figs. 5 and 6, respectively, in 
a finite range of p and 1/p. At p = 1, perfect agreement 
with the MC data by Mon [U [full circle in Figs. 5(a) 
and 6(a)] is found. 

Figs. 5 and 6 demonstrate the weakness of the n de- 
pendence at T c . On the basis of the monotonicity of the 
curves for the case n — oo (thin curves in Figs. 5 and 
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FIG. 8: (Color online) Critical amplitudes (a) $(0,p), (p39]l . 
and (b) E(0, p), (|5.11|) . at T c in three dimensions as a function 
of the inverse aspect ratio 1/p for n — 1 [thick lines, from 
(15. ip ] and in the large-n limit [thin lines, from (|3.6p ]. The 
maximum —0.3658 of the n = 1 line in (a) is at (l/p) ma x = 
0.3223. The dashed lines are the extrapolations of the n = 1 
lines from (l/p) m ax to 1/p = corresponding to cylinder 
geometry. The dotted lines represent (|5.1|> m the regime of 
small 1/p < (l/p)max where our perturbation theory is not 
applicable. See also Fig. 6. 



in Figs. 7(a) and 8(a). This leads to our prediction of the 
n = 1 amplitudes of the scaling functions of the excess 
free energy density at bulk T c for the film and for the 
cylinder in three dimensions: 

F fUm (0) « F(0, p = 1/4) = -0.164, (5.8) 
$ c „*(0)«$(0,p = 3) = -0.366. (5.9) 

The corresponding results for the Casimir amplitudes 
X(0,p) and S(0,p) are shown in Figs. 7 (b) and 8 (b); 
they follow from those of F(Q,p) and &(0,p) by means 
of the exact relations [compare (|2.29l) and (I2.41[) ] 



X(0,p) = (d-l)F(0,p)-p 



d F(0,p) 
dp ■ 



3(0, p) = -*(0,p) + (1/p) 



d(l/p) ■ 



(5.10) 



(5.11) 



From (|5.8p and (|5.9p we obtain our prediction of the 
n = 1 amplitudes of the Casimir force scaling functions 
at bulk T c for the film and for the cylinder in three di- 
mensions [dashed lines in Figs. 7 (b) and 8 (b)]: 

Xf ilm (0) = X(Q, 0) = 2F film (Q) = -0.328, (5.12) 
E cy i(0) ee 5(0, oo) = -$ cy i(Q) = 0.366. (5.13) 

Our results for Ffu m (0) and Xfu m (0) are in good agree- 
ment with the MC estimates [lj| Ap = —0.152 and 
2A P = -0.304 [triangles in Fig. 7] for the three- 
dimensional Ising model in film geometry at bulk T c . The 
previous e expansion results up to 0(e) [19] [squares in 
Fig. 7], and up to 0(e 3 / 2 ) [2(| [diamonds in Fig. 7], are 
in less good agreement with the MC estimates. 

It would be interesting to test our predictions for $ cy ;(0), 
(|53|) . and E cy i(0), (|5TT5|) . by MC simulations for the 
three-dimensional Ising model in cylinder geometry. 



6) we expect monotonicity also for the n = 1 curves. As 
shown in the magnified plots of Figs. 7 (a) and 8 (a), 
F(0,p) and Q(0,p) indeed have the expected monotonic 
behavior, but only in the restricted range p > p m ax = 
0.2470 and 1/p > (l/p) max = 0.3223, respectively. As 
expected on general grounds, the lowest-mode separa- 
tion approach should fail for sufficiently small p < p m ax 
or 1/p < (l/p) m ax, respectively, near the film and the 
cylinder limit (dotted lines in Figs. 7 and 8) where the 
higher modes are no longer well separated from the single 
lowest mode. Thus our hypothesis of monotonicity pro- 
vides the following quantitative estimate for the range 
of the aspect ratio p within which our lowest-mode sep- 
aration approach for the free energy is expected to be 
reliable: 

I/ 4 ^P^ 3. (5.7) 

Furthermore we expect a negligible dependence on p and 
1/p in the range p < p max or 1/p < (l/p) max , respec- 
tively, corresponding to the extrapolations (dashed lines) 



B. Finite-size scaling functions 

Now we turn to a discussion of the temperature de- 
pendence. In Figs. 9 and 10 we show the scaling 
functions F ex (x,p), X(x,p), & ex (x^,p), and 3(xu,p) for 
ii = 1 in three dimensions for slab, cube, and rod geome- 
tries, respectively, with finite aspect ratios in the range 
1/4 < p < 5/2 , as derived from (14351) . QZTQ , (l2~^9l . 
(j2.40p , and (]2.4ip . It is expected that these curves are 
applicable to the central finite-size regime |x| < 0(1) 
and |ajii | < 0(1) but not to \x\ ^> 1 and | 3> 1. (For a 
more precise estimate see below.) Figs. 9 and 10 should 
be compared with the corresponding Figs. 3 and 4 for 
the case n = oo. 

We see that there are significant differences between the 
cases n — 1 and n — oo. Figs. 9 (a) and 10(a) exhibit 
a nonmonotonicity of F ex (x,p) and & ex (x\\,p) for n = 1 
with minima slightly below T c for all p. Such minima 
should also persist in the n = 1 film (p = 0) system and 
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FIG. 9: (Color online) Scaling function (a) F ex (x,p), (2~27)l , 
(035), and (b) X{x,p), (2ZD, (J227), (£55) , as a function 
of x = i(Z,/£ +) 1/ " for n = 1 in three dimensions for slab 
geometry with finite aspect ratio p = 1/4 (solid lines), p — 1/3 
(dotted lines), p — 1/2 (dashed lines), p = 2/3 (dot-dashed 
line), p = 1 (double-dot-dashed line). Thin lines in (a): e 
expansion results for p = from [Tg| . [2lj . 



shown by the double-dot-dashed curves in Figs. 9 (b) and 
10 (b) [see also Figs. 5 and 6]. In addition, the Casimir 
force for n = 1 in a cube changes sign at x — —0.884 
and is negative for x < —0.884, contrary to the case 
n = oo below T c [Figs. 3 (b) and 4 (b) for p = 1]. Thus 
our theory predicts that, in a cube, there is only a small 
positive region between x = —0.884 and 5 = 0. 



Or 




-0.4- - 



-0.6 



-0.8 



l/p= 2/5 

1/p = 1/2 

1/p = 2/3 

1/p = 1 




in the n = 1 cylinder (1/p — 0) system whose scaling 
functions should be close to our curves for p = 1/4 and 
1/p = 2/5, respectively. There is no good agreement at 
T c between our p = 1/4 curve in Fig. 9 (a) and the e 
expansion results (thin lines ) of [l9l. r2lT| for p = 0. The 
latter exhibit an unphysical singularity at x = (i.e.. at 
bulk T c ) that arises from the e expansion results [19|, 
for the term F(x, p) in (|2.27l) which should be an analytic 
function of x near x = since the film transition occurs at 
a distinct temperature T c ju m below bulk T c . Our curves 
contain a different type of singularity at x = that arises 
from subtracting the singular bulk part F^(x) in (|2.27l) ; 
this singularity is very weak and not visible in Figs. 9 
and 10. 

In Fig. 9 (b) our results show an unexpected structure of 
the Casimir force scaling function X near bulk T c where 
local maxima occur with increasing p > 1/4. The small 
shoulder for p = 1/4 was already noticed previously [24| . 
This structure with local maxima does not exist for n — 
oo. Such maxima also persist in the regime of p > 1 as 
shown in Fig. 10 (b). As a special feature of the case 
p = 1, X and 5 vanish at bulk T c in three dimensions, as 



FIG. 10: (Color online) (a) Scaling function (a) & ex (x\\ , p), 
(230), (12371) . (435), and (b) S(x\\,p), (2^41) . as a function 
of Sii = t(Z/||/£o+) 1 ^ ! ' for 7i = 1 in three dimensions for rod 
geometry with finite aspect ratio 1/p = 1 (double-dot-dashed 
lines), 1/p = 2/3 (dotted lines), 1/p = 1/2 (dashed lines), 
1/p = 2/5 (solid lines). 

On purely theoretical grounds, it is difficult to provide a 
precise quantitative estimate for the range of validity of 
our perturbation approach with regard to the dependence 
on the scaling variable x. Valuable information, however, 
has been made available to us by Hasenbusch (36| who 
performed MC simulations for the free energy density of 
the three-dimensional Ising model in a cubic geometry. 
These data are shown in Fig. 11, together with our theo- 
retical curve derived from Q2.27P and (|4.55[) . We see that 
there is good agreement in the range —0.05 ^5^3 but 
significant deviations exist well below T c ; small but sys- 
tematic deviations exist also well above T c . In particular, 
our perturbation result for F ex has an algebraic approach 
to a finite limit F ex (oo,p) for x — > oo whereas there 
should be an exponential decay towards zero (see Sec. 
VI). From this comparison it is obvious that the lowest- 
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mode separation approach needs to be complemented by 
a perturbation approach that is valid outside the central 
finite-size regime. Such an approach will be presented in 
the subsequent section. 

Additional valuable information comes from a compari- 
son of our Casimir force scaling function with earlier MC 
data for periodic b.c. in the small - p regime [lj|. We re- 
call that the lower limit of applicability of our calculation 
is p = 1/4 and that the Casimir forces should depend only 
weakly on p for p < 1/4, thus it is reasonable to com par e 
our result for p = 1/4 with MC data for p = 1/6 [16j . 
This comparison is shown in Fig. 12 . Also shown are 
the previous e expansion results for p = from [l|| [2l[ 
which exhibit the same kind of singularity at x = as 
in Fig. 9 (a). We see good agreement of the MC data 
with our fixed d perturbation theory in the whole range 
— 2 < x < 20. There are systematic deviations only for 
x < — 2 which are less pronounced than those for F ex 
in the same region. In the subsequent section we shall 
explain this different degree of agreement between our 
theory and the MC data shown in Figs. 11 and 12. 
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FIG. 12: (Color online) Scaling function X(x,p), (229]), 
HUT]) , (035), for n = 1 as a function of x = t(L/£ +) 1/ " 
in three dimensions for slab geometry with p = 1/4 (solid 
line) and MC data for the d = 3 Ising model with p — 1/6 by 
Vasilyev et al. [l(|. Thin lines: e expansion results for p = 
from [HHH]. See also Fig. 14 (a). 




FIG. 11: (Color online) Scaling function F BX {x,l), (223), 
(|4.55p . for n = 1 as a function of x = t(L/£o+) 1 ^" in three 
dimensions for cubic geometry (solid l ine) and MC data for 
the d = 3 Ising model by Hasenbusch [3g|. See also Fig. 13 
(c). 



VI. Perturbation theory outside the central 
finite-size regime 



Outside the central finite-size regime, there is no need 
for separating the lowest mode, thus ordinary perturba- 
tion theory with respect to uq should be appropriate. 
By "outside the central finite-size regime" we mean the 
regions below the dashed lines in Fig. 2. Within these re- 
gions, it is necessary to further distinguish between scal- 
ing and nonscaling regions (the nonscaling regions cor- 
respond to the shaded regions in Fig. 2, see also Fig.l 
of Q). The central parts of both the scaling and the 
nonscaling regions still belong to the asymptotic critical 
region \t\ <C 1 and L s> a, L\\ ^> a. 

Here we perform the corresponding analysis at the one- 
loop level. In subsections A - C, we shall consider 
the scaling region outside the central finite-size regime. 
The total scaling region can be roughly characterized 
by L/a > 1, L\\/a > 1, f±/5 > 1, and L/£± < 
24(£±/5) 2 ,£||/£± ^ 24(£ ± /5) 2 where £± is the second- 
moment bulk correlation length above and below T c , re- 
spectively. (Note that this characterization also includes 
the central finite-size regime which is part of the total 
scaling region.) The latter restrictions follow from the 
conditions (16. 18)) for the nonscaling regions that will be 
studied in subsection D below. In order to distinguish the 
perturbation results of this section from those of Sees. IV 
and V we use the notation f£_ loop , f^L loop , etc. 
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A. Perturbation theory well above T c 

Ordinary perturbation theory for the excess free energy 
density (|2.9p for n = 1 above T c yields in one-loop order 

f?-'L P = ^ E Mfa - ^ + ^(k)]a 2 } 

k 

- iyin{[ro-r 0c + ^(k)]ci 2 }. (6.1) 

k 

Here we have already replaced ro by r$ — r$ c which is 
justified since r 0c ~ O(u ) [see (|4.28l) ]. Because of the 
k = term, the sum exists only for ro — r$ c > 0. The 
evaluation of the excess free energy density is outlined in 
App. A. 

In the scaling region in 2 < d < 4 dimensions, the large- 
k dependence of SK(k) does not matter and the leading 
contribution is obtained by taking the continuum limit 
a — >• at fixed ro — ro c > 0. For a renormalization-group 
(RG) treatment in the scaling region see (10.5) - (10.13) 
of Q. Neglecting nonasymptotic corrections to scaling 
we obtain the scaling function 

F?-\ + oo P (x,P) = \ Go(x 2u ,p) + 0(u*) , (6.2) 

where Qo is given by (|3.7p and x" = £+ = £o+£~ !/ 

(for £ + see (5.16) of Q)- For large x, F^\~^ op decays ex- 
ponentially to zero according to the asymptotic behavior 

V(rf-l)(^) V exp(-ip , (6.3) 

apart from corrections of 0(e~ 2x , e~ 2a: M), with xY, = 
x v jp. Eq. ((OJ) follows from (|A~T6)l in App. A for 
£,+ /a 3> 1. We see that the scaling variable Xn appears 



in a natural way in the second term of (|6.3[) . For p = 1, 
([EgD and agree with Eqs. (10.10) - (10.12) of 0. 

The corresponding scaling functions ^^2io 0p -, X^_ loop , 

^j-loo p and ^Zyti p, X ^s ymp, ^tsymp follow from ([621), 

(jO]) . IH2S]), ([2T4TJ|) . and (|2T4Tj) . respectively. 



B. Perturbation theory well below T c 

Perturbation theory for bulk quantities below T c within 
the ip model for n = 1 at vanishing external field h may 
be formulated by first starting with the perturbation ex- 
pression at finite external field h > (or h < 0) and at 
finite volume V — L^ 1 L, then performing the thermo- 
dynamic limit V — > oo at finite h > (or h < 0), and 
subsequently performing the zero-field limit h — > 0+ (or 
h — > 0_). Applying this procedure to the free energy 
density f(t,L,L\\,h) implies that only the contributions 
of a single bulk phase with a positive (or negative) spon- 
taneous bulk magnetization are taken into account in the 
calculation of 

fb(t) — lim lim f(t,L,Lu,h)— lim lim f(t,L,Lu,h 

In MC simulations of finite Ising models at vanishing ex- 
ternal field, however, all configurations of both phases 
with positive and negative magnetization do contribute. 
In this case, the order-parameter distribution function 
has two finite peaks with equal heights in the positive 
and negative ranges of the magnetization [171, |25j, |4a, |49( . 
For T — > 0, these two peaks are well separated. In or- 
der to account for this fact in an analytic treatment of 
the model well below T c , it is appropriate to formu- 
late perturbation theory such that an expansion is made 
around the two separate peaks of the order-parameter 
distribution function that exist at h = 0. 



In the following we perform this approach at the one- 
loop level in order to calculate fl^^, op {t, L, L|| , h = 0) 
well below T c . First we decompose the lattice variable 
ipi of the Hamiltonian H, (|2.1j) . as ipi = M > , m f + 5ipi 
with the positive mean-field order parameter M >>m f — 
[— r /(4u )] 1//2 > 0. Keeping only the Gaussian terms 
of H up to 0[(Sifi) 2 ] corresponding to a one-loop ap- 
proximation leads to the dimensionless partition function 



(compare (Bl) of 0) 

Z> = exp \V-$-\ TT ( ^ ] ■ (6.5) 

16 Uo J1 k L \ v [-2ro + ^(k)]5 2 y 

It can be rewritten as 

Z> = exp(-I// fc -) Z~ (6.6) 

where fr is the bare bulk free energy density (14. 18|) in 
one-loop order below T c and 

ZSf=exp{V7 ea; } (6.7) 



19 




FIG. 13: (Color online) Scaling functions F ex (x,p), (a) - (c), and Q ex (icy , p), (d), as a function of x = t(L/^o+) 1 ^ 1 ' and 
xn = t(L||/^o+) 1// ", respectively, for n = 1 in three dimensions for several values of the aspect ratio p. Thick lines: improved 
perturbation theory according to (|4.55[) and (|2.40[) . Thin lines: one- loop perturbation theory according to (16. 2[1 and (|6.10|) . 
The thin lines diverge for t — > 0. The asymptotic value of the thin lines is — p 2 In 2 for x — s> — oo in (a) - (c) and —(1/2) In 2 for 
a; II — > — oo in (d). MC data in (c) for the d — 3 Ising model by Hasenbusch [361 ]. 



is the finite-size part of Z > with the contribution 
f BX = ^7E ln {[-2r + ^(k)]5 2 } 

- ^J\n{[-2r + 5K(k)]a 2 } (6.8) 

k 

to the excess free energy density below T c . The partition 
function Z>, however, is incomplete with regard to the 
finite-size contributions as it does not take into account 
the fluctuations around the negative mean-field order pa- 
rameter M <tin f = — [— ro/(4uo)] 1 / 2 < 0. A decomposi- 
tion of ifi as ifi = M <im f + 6tfi and an expansion of H up 
to 0[(d<fi) 2 ] leads to the one-loop partition function Z < 
which is of course the same as Z > . Thus the total finite- 
size part Zf^ + Z^ of the partition function in one-loop 
order is then given by 2 exp{V f ex }. The corresponding 
total excess free energy density is 

A^W = + ^E ln {h2(r -r 0c ) + ^(k)]a 2 } 

k 

^J\n{[~2(r -r 0c ) + SK(k)}a 2 }. (6.9) 

k 



Here we have again replaced — 2tq by — 2(ro — ro c ) in the 
spirit of perturbation theory up to 0(1). The result (|6.9[) 
is identical in form with the (corrected 32] ) result derived 
previously for cubic geometry Q- The non-exponential 
finite-size term — T/ _1 ln2 is known from previous work 
on finite-size effects in Ising models in a block geometry 
of volume V [25j . Thus this term is not specific for the 
n = 1 (p 4 theory but rather general for systems with a 
two-fold degeneracy of the ground state. According to 
the definition of the Casimir force (|2.10[) . the constant 
term —(In 2)/V in (|6.9[) does not contribute to Fcasimir- 

The derivation presented above is, of course, not exact 
but is valid only well below T c where the two peaks of 
the order-parameter distribution function at h = are 
well separated and where their wings do not overlap sig- 
nificantly. 

The evaluation of (|6.9p as well as the RG treatment 
are parallel to that for fi^ op above T c . Neglecting 
nonasymptotic corrections to scaling we obtain the scal- 
ing function in the scaling region well below T c 



nXov&P) = - ^ 2 + - e (|2i| 2l » (6.10) 
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FIG. 14: (Color online) Scaling functions X(x,p), (a) - (c), and n(x«,p) t (d), as a function of x — t(L/^o+) 1 ^ v and x\ 
t(I/||/fo+) 1 ^" 5 respectively, for n = 1 in three dimensions for several values of the aspect ratio p. Thick lines: improved 
perturbation theory according to (|4,55[) , (|2,29|l , and (12.41 [) . Thin lines: one- loop perturbation theory according to (I6.2|l and 
(16.10[) . The thin lines diverge for t — ¥ 0. MC data in (a) for the d = 3 Ising model with p = 1/6 by Vasilyev et al. [16 



with L/£- = \2x\ v where 



.11) 



is the bulk second-moment correlation length below T c 
and where Go is given by (|3 . T[) . In contrast to the van- 



ishing of i^'foop for x — > oo, F^_' loop approaches a finite 
value — /9 d ~ 1 ln2 for x — > — oo, as noted already in [32j . 
According to (|6.10l) and (|6.3[) . this approach has an ex- 
ponential form described by the asymptotic behavior 



F ex- 
asymp 



-A 1 In 2 



V 2vr 



12X1,1^^ 



- exp(-|2£,||' ) 

apart from corrections of 0(e~ 2 ' 2:E ' 



exp(-|25| I/ ) 
(6.12) 

")■ 



X ~ 

loop' 1 — loopi 



The corresponding scaling functions $]_' 

^ ^ ^ Smf 4"s;, mp. B qsymp foll ° W fr ° m <l 6 - 10 l>» 

([5TT^]) . (j2~29| . (pOO]) . and (pTITj) , respectively. 



As noted above, the constant term — p In 2 does not 
contribute to X. This explains why the perturbation re- 
sult for X of Sec. V (as shown in Fig. 12) is in better 
agreement with the MC data below T c than the corre- 



sponding result for F ex shown in Fig. 11. Thus, in con- 
trast to F ex, ~ which has a finite low-temperature limit 
— p d_1 ln2, our theory predicts that the Casimir force 
scaling function has an exponential decay towards zero 
for i — > — oo. From (|6. 12|) and (|2.29[) we obtain 



(tl + |2ir) (ffi^ eip( _ |2ir) 



+ P d (d-l) 



\2x\\ 



2n 



exp(-|2i||r). (6.13) 



It is suggestive to expect that the formulae (|6.3|) , (|6.12[) , 
and (|6.13p are applicable even to d = 2 - dimensional 
systems. It would be interesting to check this point for 
the example of the two-dimensional Ising model with pe- 
riodic b.c. in rectangular geometry. 



C. Predictions for the whole scaling region 



On the basis of the three perturbation results (|4.55[) . 
(|6.2[) , and (|6.10p we are now in the position to present 



21 



quantitative predictions for the various scaling functions 
over the whole range of the scaling variables —15 < x < 
20 and —15 < X\\ < 20. These scaling functions are 
shown in Figs. 13 and 14 for various values of the aspect 
ratio p in three dimensions. 

The thin lines are based on one-loop perturbation the- 
ory (|6.2p and (|6 . 10[) and are applicable only away from 
T c outside the central finite-size regime. For T — > T c , 
one-loop perturbation theory breaks down which implies 
that the thin lines diverge for t — > 0. The thick lines are 
based on our lowest-mode separation approach presented 
in Sees. IV and V which is applicable to the central finite- 
size regime including T — T c . This improved perturba- 
tion approach provides a bridge through T — T c between 
the simple finite-size critical behavior represented by the 
thin lines well away from T c . The lowest-mode separa- 
tion approach is not applicable, however, to the regions 
|x| ^> 1 and |x||| 3> 1. Our Figs. 13 and 14 demonstrate 
that one-loop perturbation theory and improved pertur- 
bation theory complement each other and match reason- 
ably well at intermediate values of the scaling variables. 
No perfect matching can be expected because of missing 
0(u*) terms in the one-loop results. Comparison with 
the MC data in Figs. 13 and 14 shows that the improve- 
ment achieved by the one-loop results is clearly visible 
in the range x > 4 and x < — 1 [in Fig. 13 (c)] and in 
the range x < —2 [in Fig. 14 (a)]. On the whole, we 
consider the good agreement of our theory with the MC 
data over the entire scaling regime — 15<5<20asa 
major success of our strategy employing three different 
perturbation approaches. Comparison with MC data for 
other values of p would be interesting. 



with the nonuniversal function 



Ak{i±,L,L\\,a) 



1 

1? 



d-1 



a 
2£I 



2tt£± 
( L \\ 



exp 



L 



exp 



(6.16) 



where 



arsinh 
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3.17) 



are the exponential ("true") bulk correlation lengths 0, 
IH, H|J above (+) and below (-) T c , respectively. This 
result applies to the regions well below the dashed lines 
in Fig. 2 including the shaded regions. Note that no 
condition is imposed on the value of < a/£± < oo 
other than that L/£± and Lu/£± are large. For L = Lu, 
(16.161) reduces to the previous result for cubic geometry 
[TJ, l50j . As a nontrivial relation between bulk properties 
and finite-size effects [33j], the lengths £ e ± describe the 
exponential part of the bulk order-parameter correlation 
function [34| in the large-distance limit in the direction 
of one of the cubic axes at arbitrary fixed T ^T c above 
and below T c (for n = 1), respectively. This relation is 
exact in the large-rt limit above T c 



It has been shown 0,113 that; because of the exponential 
structure of the function A±, the a dependence of £ e ± 
cannot be neglected even for small 5/£± <C 1 if 



L>24&/a 2 , i||>24^/a 2 



(6.18) 



D. Exponential nonscaling region 



are sufficiently large. The conditions (|6. 18|) follow from 
the second term in the expansion of the function (|6.17[) 
for small a/£± 



So far we have eliminated the dependence on the lattice 
spacing 5, by taking the continuum limit. In earlier work 
it was pointed out for confined systems in an L d geometry 
0, |4(| and in film geometry Q that the finite lattice con- 
stant a becomes non- negligible in the limit of large L/a 
at fixed T ^ T c in the regime where the finite-size scaling 
function has an exponential form. The same arguments 
apply to the present system in a finite block geometry. 
As shown in App. A, the excess free energy density in 
one-loop order attains the following form in the limit of 
large L/a, large Ln/a, large L/£±, and large £||/£± 



tex,+ —A(C 
J asymp \Ss-i 



,L,L||,o), 



(6.14) 



!■( .1 

J as 



In 2 . ,„ T T 
-— +-A-(£-,L,Li\,a), 



(6.15) 



-f-V 

24 Ve±/ 



(6.19) 



appearing in the exponenti al p arts of the function 
•A±(£±iL,Lu,a) (see also [33|,|46[). The second term in 
(|6.19p is not negligible even for small a/£± "C 1 if the con- 
ditions ()6.18|) are satisfied. This implies that finite-size 
scaling and universality are violated in the large |i| and 
tail of L d f^ x at any o/£± > even arbitrarily close 
to T c because ultimately, for |5| —t oo and \xu \ — > oo (i.e., 
for large L and Ln at fixed \t\ > 0), the tail of L d f^ x be- 
comes explicitly dependent on a. As shown in Sec. X of 
, the tail depends even on the bare four-point coupling 
Mo through £±: strictly speaking it is even necessary to 
keep the complete nonasymptotic (uq dependent) form 
of £± at finite a. Thus no a - independent finite-scaling 
form with a single scaling argument cx tL x l v can be de- 
fined in this exponential large \x\ and large region. 
Higher-loop contributions cannot remedy this violation. 
The same reservations apply, of course, to the critical 
Casimir force and its scaling form. 
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Note added: The predictions of Ref. [32j and of the 
present paper are in good agreement with recent Monte 
Carlo data for the three-dimensional Ising model by A. 
Hucht, D. Griineberg, and F.M. Schmidt, Phys. Rev. E 
83, 051101 (2011). 
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Appendix A: Gaussian free energy 



We consider the Gaussian model, i.e., the Hamiltonian 
(|2.ip for uq = 0, and calculate the excess free energy 
density in a rectangularLjj -1 x L geometry. This calcu- 
lation will lead to the evaluation of the sums in (13.11) and 
([372]) as well as to the derivation of ((47231) - ([4727) , ([672]) . 
([6731) . ([BTTO]) . ([6712]) . and (foTTB]) . Since the calculation is 
largely parallel to that of [7| we skip some of the details 
of the derivation. 



The Gaussian excess free energy density per component 
divided by k B T is f^ auss = ±A(r ,L^L, K iyJ ,a), 

A(r , L|| , L, Kij,a) = V' 1 ^ Mfa> + 6K(k)]a- 2 } 

k 

ln{[r Q + SK(k)]~a 2 } (A.l) 



where the sum J^k an d the integral J k have finite cutoffs 
±n/a for each k a . Using the Poisson identity [33|, [47j we 
obtain the exact representation 



A(ro,L\\,L,Kij,a) = - dyy 1 



-r a y 



exp{— SK(q,p)a 2 y + iq ■ min + ipnL}(A.2) 



with q-m = 5TJ a =i a a rn a where Yl' m n means summation 
over all integers m = (mi,m2, ...,md-i) and n without 
the single term with m = 0,n = 0. In the following we 
evaluate A for L 3> a, and Lu 3> a in two regimes. 



1. Central finite-size regime 



l ll 

We assume large L/a, large iy/ffl, small < r a -C 1 

and fixed < Lr]J 2 < 0(1), < L^ 2 < 0(1) which we 
refer to as the central finite-size regime. In this regime, 



the large - (q,p) dependence of 5K(q,p) does not mat- 
ter. Therefore we may replace 5K(q,p) by its long - 
wavelength form (I2.20|) and let the integration limits of 
/ and of J go to oo. This leads to the scaling form of 
the Gaussian excess free energy 

f Gauss = \^r Q ,L h L) = h- d G (r L 2 ,p) (A.3) 



where Qj(r L 2 ,p) is defined in (|3.7I) . Interpreting (|A.3[) 
as a one-loop contribution of the ip 4 model and applying 
the renormalization procedure parallel to that described 
in Sec. X A. of 0| we arrive at the one-loop scaling 
function presented in 



The function Go(roL 2 , p) diverges for tqL 2 — > which 
comes from the large- z behavior of K(z) « 1 in the last 



term 



We find 



pK(p 2 z) K(z) « p d ~ x of the integrand of (j3"7Tl) . 



Go(r L 2 , p) « p*" 1 In (DL_) + C (p) (A.4) 

for roi 2 <C 1. In order to determine the con- 
stant Cq(p) we add and subtract the divergent term 
p' 1 ^ 1 ln[(roL 2 /(47r 2 )] by rewriting Qo(roL 2 , p) in the form 



g (r L 2 ,p) = p d - 1 ln( 



r L 2 
V 4tt 2 



dz 

z 



exp 



r L 2 z 
4tt 2 



7T\ d /2 
2 / 



pK(p 2 z)] d 1 K{z) + p d - 1 



p^e-" 



(A.5) 

The integral in (IA.5[) has a finite limit for roL 2 — > which 
yields the constant 



Co GO 



d-1 



(- 



d/2 



- P K{p 2 z)\ K(z) + p d -\\-e- z ) 

Eq. (IA.4I) implies that the function 

g 1 (r L 2 1 p) = 
has the divergent behavior 



dg (roL 2 ,p) 



Gx{r Q L 2 ,p)n- 



P 



r L 2 



(A.6) 



(A.7) 



(A.8) 



for r L 2 < 1. The asymptotic behavior (|A.4[) and (|A.8[) 
is needed in the discussion of the low-temperature limit 
in Sec. III.C . 

The function Go(roL 2 , p) decays exponentially for large 
tqL 2 3> 1. From (|A.16[) we obtain for rod, 2 <C 1 and 
r L 2 > 1 

- 2X (d-l)/4 



g {r L 2 ,p) » -2 



4tt 2 

2 \ (d-l)/4 



exp(-Lro /2 ) 



exp(-iry 2 /p)- (A.9) 
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For p=l, Eq. (|A~9)) agrees with Eq. (10.12) of 



The result (|A.3|) is sufficient to derive the higher-mode 
sums Si(r L,L 7 p), f|4. 1 1[) and f|4. 12[) in the central finite- 
size regime. We obtain 



i^ln{[ro + ^(k)]a 2 } 



ln{[r -f-«f(k)]5 2 } 



2 4tt 2 



+ 1 £ P d ln(r a 2 )+ l^J ( ro L 2 ,p) (A.10) 

with Jo (x 2 > p) defined by (|4.26l) . By means of differentia- 
tion with respect to ro we obtain from (|A.10[) for m = 1, 2 

v ^ i>° + ^( k )] _m = /Vo + ^(k)]- m 



d-1 



-(*>)" 



2m — d 



-I m (r L 2 ,p) (A.ll) 



(4tt 2 )™' 

with I m (roL 2 , p) defined by l|4.27p . For the bulk integrals 
see 0- 



2. Exponential regime above T c 



1/2 ' 

Now we assume LrJ ^> 1 and Lur^'^ ^> 1 at finite 



l li 

p = L/L\\ for fixed r a > which we refer to as the 
exponential regime since A(ro, L\\ , L, Kij, a) will attain 
an exponential L and L|| dependence in this regime. In 
this regime the complete k dependence of the microscopic 
interaction 8K(k) does matter. We use the nearest- 
neighbor interaction (|2.19l) in the form 

2 2 
SK(q,p) = — [! - c °s(ag Q )] + — [1 - cos(ap){A.12) 

a— 1 

Generalizing the derivation of to block geometry we 
obtain for large L/a and L\\/a but arbitrary ro = roa 2 > 

(compatible with Lr\j 2 3> 1 and L\\r\j 2 ;§> 1) 



II 'o 



1/2 



dz 

a d (2nL/a) d / 2 



2{d- 1) 



cxp[$(z,f )£/q] 

z (d+l)/2 g l/2 

exp[$(z,f )£||/a] 



a d (2^L\s/~a) d l 2 J z (rf+l)/2 g l/2 




(A.13) 



where q = (1 + z 2 ) 1 / 2 with 

$(z, f ) - -(1 + r /2)z + g + In ) . (A.14) 

VI + <?/ 

The maximum of the function $(z, fo) in the exponential 
parts of the integrand of (IA.13[) is at z — z where 



ro 



-1/2 



(A.15) 



Expanding $(z, fo) around z = z up to 0[(z — z) 2 ] and 
performing the integration over z we finally obtain the 
Gaussian excess free energy density for large L/a and 
large Ln/a at arbitrary fixed ro > 



rex 

I Gauss ,asymp 



(d-l)/2 



(d-1) 



L || 



:x (d-l)/2 



2ttz 



(A.lf 



with the exponential ("true") bulk correlation length of 
the Gaussian model 



We recall that r, 



a 
2 

-1/2 



arsinh [ — 



(A.17) 



qi is the second-moment bulk 
correlation length of the Gaussian model above T c . For 
L = L|| (cube), (|A.16|) yields the previous result of Eq. 
(B24) of 0. No universal finite-size scaling function of 
the Gaussian model can be defined in the region L > 
24(£f) 3 /a 2 and £j > 24($g) 3 /a 2 because of the explicit 



a dependence of 



(IAT61) 



and ([A~T7| . 



Within a RG treatment of the ip 4 lattice model the Gaus- 
sian results (fOt and (|A~T6)) can be considered as the 
bare one- loop contributions to the excess free energy den- 
sity. By means of such a RG treatment at finite lattice 
constant a parallel to Sect. 2 and App. A of (46|, these 
results acquire the correct critical exponents of the n = 1 
universality class including corrections to scaling. This 
leads to the one-loop results at finite a in Sec. VI D. 
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